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Abstract. For N interacting particles in a one dimensional random potential, we study the structure of 
the corresponding network in Hilbert space. The states without interaction play the role of the "sites". 
The hopping terms are induced by the interaction. When the one body states are localized, we numerically 
find that the set of directly connected "sites" is multifractal. For the case of two interacting particles, 
the fractal dimension associated to the second moment of the hopping term is shown to characterize the 
Golden rule decay of the non interacting states and the enhancement factor of the localization length. 
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The wave functions of one particle in a random po- 
tential have been extensively studied. In two dimensions 
within the localization domains [1] the large fluctuations 
of their amplitudes have a multifractal character. In one 
dimension, the elastic mean free path I and the localiza- 
tion length L\ coincide, preventing a single one particle 
wave function to be multifractal over a significant range 
of scales. The description of the correlations existing be- 
tween the localized eigenstates is more difficult. This is 
quite unfortunate, since a local two-body interaction re- 
organizes the non interacting electron gas in a way which 
depends on the spatial overlap of (four) different one par- 
ticle states. When one writes the TV-body Hamiltonian in 
the basis built out from the one particle states (eigenba- 
sis without interaction), this overlap determines the in- 
teraction matrix elements, i.e. the hopping terms of the 
corresponding network in Hilbert space. In this work, we 
numerically study the distribution of the hopping terms in 
one dimension, when the one body states are localized. It 
has been observed [2] that this distribution is broad and 
non Gaussian. We give here numerical evidence that this 
distribution is multifractal. Moreover, since the obtained 
Renyi dimensions do not depend on L\, simple power laws 
describe how the moments scale with the characteristic 
length L\ of the one body problem. Since the main ap- 
plications we consider (Golden rule decay of the non in- 
teracting states, enhancement factor of the localization 
length for two interacting particles) depend on the square 
of the hopping terms, we are mainly interested by the 
scaling of the second moment. For a size L w L\, we show 
that, contrary to previous assumptions, the iV-body eigen- 
states without interaction directly coupled by the square 
of the hopping terms have not a density of the order of 



the two-body density ^2(^1) cx L\, but a smaller density 
pf (Li) oc L ^ a ("^\ The dimension f{a(q = 2)) (« 1.75 
for hopping terms involving four different one body states) 
characterizes the fractal set of TV-body eigenstates with- 
out interaction which are directly coupled by the square 
of the hopping terms. 

We consider N electrons described by an Hamiltonian 
including the kinetic energy and a random potential, plus 
a two-body interaction: 

H = t a d+ a d aa + U ^2 Ql^d+^id^ds^ (1) 

The operators d+ a (d aa ) create (destroy) an electron in 
a one body eigenstate \a > of spin a. Noting ^ a (n) the 
amplitude on site n of the state \a > with energy e a , the 
interaction matrix elements are proportional to the QTg 
given by: 



(2) 



This comes from the assumption that the interaction 
c n| c «l c nT c ™T is local - Tne c t,] ( c «^) crea t e (destroy) 
an electron on the site n and = ^2 n ^ a {n)c^ a . When 
U = 0, the Hamiltonian is diagonal in the basis built 
out from the one particle states, and the N body states 
dliLi dt iC n l u >) can De thought as the "sites" with en- 
er gy Eili of a certain network which is not defined 
in the real space, but in the iV-body Hilbert space. When 
U 7^ 0, different "sites" can be directly connected by off- 
diagonal interaction matrix elements. Therefore, one can 
map [3] this complex A^-body problem onto an Anderson 
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localization problem denned on a particular network in the 
TV-body Hilbert space. Since the interaction is two body, 
only the "sites" differing by two quantum numbers can be 
directly coupled. This restriction will not matter [4] for 
N = 2 and (under certain approximations) may yield a 
Cayley tree topology [3] for the resulting network, if TV is 
large. We study the additional restrictions coming from 
one body dynamics. 

We summarize a few evaluations of the second moment 
(q = 2) of Cpjp which have been previously used. Case (i): 
The one body Hamiltonian is described by random matrix 
theory (RMT). The statistical invariance under orthogo- 
nal transformations 0(M) implies that ((Q^) 2 ) ~ 1/M 3 
where M is the number of one body states. Case (ii): 
The system is a disordered conductor of conductance 
g. An estimate [3] based on perturbation theory gives 
{{UQ 1 ^) 2 ) oc (A/g) 2 . Since the one particle mean level 
spacing A oc 1/M, this perturbative result coincides with 
the previous RMT results if one takes M = g 2 . Moreover, 
it is valid only if all the one particle states appearing in 
Eq.(2) are taken from a sequence of g consecutive levels 
in energy. Otherwise, Qf^a can be neglected. Case (hi): 
The system is a disordered insulator. Shepelyansky [4] 
in his first study of the two interacting particles (TIP), 
assumes a RMT behavior for the M = Lf components 
of the wave function inside the localization domain, and 
neglects the exponentially small components outside this 
domain. When the dimension d = 1, one gets a term 
(Ql^) 2 ~ l/Lf for the terms coupling a TIP state \a(3 > 
to L\ TIP states \<y& >. This estimate for g < 1 differs 
from the one valid when g > 1 under two important as- 
pects: not only M Lf instead of g 2 , but the condition 
for a large hopping is entirely different. In the insulator, a 
large hopping term is not given by four one particle states 
close in energy, but by four states close in real space, i.e. 
located inside the same localization domain. Ponomarev 
and Silvestrov have criticized [5] this estimate, using an 
approximate description of a localized state for weak dis- 
order. They note that the density of TIP states coupled 
by the interaction is sensibly smaller. 

For a more accurate study of Q^a in one dimension, we 
consider a spin independent one particle Anderson tight 
binding model with L sites and nearest neighbor hopping 
(t = 1). The on-site potentials V n are taken at random 
in the interval [-W,W] and the boundary conditions are 
periodic. L\ is estimated from the weak disorder formula 
L x w 25/ W 2 . The |Q^| are calculated using Eq.(2) and 
numerical diagonalization of the one particle Hamiltonian. 
QTg, for fixed a and (3 is a two-dimensional object which 
is not defined in the real 2d space, but in the space of two 
one particle quantum numbers 7 and <5. Those states I7 > 
(and I <5 >) can be ordered in different ways: (a) spectral 
ordering by increasing eigenenergy, (b) spatial ordering 
by the location n 7 of their maximum amplitude, from one 
side of the sample to the other, (c) momentum ordering if 
W = 0. Let us note that ordering (b) becomes meaningful 
only in the localized regime (L > L±). 



We first study the matrix element Q2? aaa > character- 
izing two electrons with opposite spins in the same state 
\ctQ > hopping to an arbitrary state |7<5 >. Hopping is very 
unlikely over scales larger than L\. The L\ large values of 
the hopping term are concentrated inside a square of size 
L 2 , as shown in Fig.l for a given sample using ordering 
(b) and a rainbow color code. Fig.l is not homogeneously 
colored, but exhibits a complex pattern which reminds us 
another bi-dimcnsional object: the one particle wave func- 
tion in a two dimensional disordered lattice. This suggests 
us to analyze its fluctuations as for the 2d one body states, 
and to check if this pattern is not the signature of a mul- 
tifractal structure. 
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Fig. 1. \QZ a a \ with |qo > taken in the bulk (n ao — 50) of a 
sample of size L = Li = 100. Spatial ordering (b). The color 
code goes from red (small value) to violet (large value) through 
yellow, green and blue. 

In analogy with the 2d one body problem, we do not 
expect that this multifractality will be valid in the whole 
(7,5) Hilbert space, but only in a limited but parametri- 
cally large domain. 

We proceed as usual (see references [6,7]) for the mul- 
tifractal analysis. For L\ and L fixed, we divide the plane 
(7,(5) into (L/D) 2 boxes of size D and we calculate the 
ensemble averaged function for different values of q 

I q (D) = ( £ \Ql s oao \\ . (3) 

i=1 \7,5eboXi / 

The existence of a multifractal measure defined in the 
(7, <5)-plane by the interaction matrix elements is estab- 
lished in the next figures. In Fig. 2, a single sample has 
been used and power laws I q (D) oc D T ^ are obtained 
over many orders of magnitude for different values of q. 

The limits of validity of these power laws are shown in 
Fig. 3. 
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Fig. 2. Power laws showing that Fig. 1 corresponds to a 
multifractal measure in the (7, 5) two dimensional plane. I q (D) 
are calculated for a single sample with L\ — L — 2500. ao in 
the bulk of the spectrum. The states I7) and S) are ordered by 
increasing eigenergy (ordering a). The dashed line corresponds 
to the RMT prediction (case (i)). 

On the left side, spatial ordering (b) for different val- 
ues of L\ is used for the states (7, 5). One can see that 
Iq(D) oc D T ( q \ for scales 1 < D < Li, as indicated by the 
arrows. The lower scale is given by the lattice spacing of 
the (n y ,ns) network in Hilbert space. The upper scale L\ 
is the largest scale compatible with a spatial overlap of the 
states 7 and S, for a fixed ao- This means that the multi- 
fractality of the interaction matrix elements QZ S oao m the 
two dimensional Hilbert space (7, <5) has the same para- 
metrically large range of validity as the one body wave 
function [1] in two dimensions (scale 1 < D < Li). Here, 
multifractality is valid for L\ matrix elements as multi- 
fractality is valid in the 2d one body problem for L\ sites. 

On the right side of Fig. 3 spectral ordering (a) is 
used for the same samples, giving the same power laws 
as with ordering (b), inside the corresponding energy 
range (A(Li) < DA(L) < 1) indicated by the arrows. 
A(x) oc x -1 is the level spacing of a segment of size x, 
and 1 is the band width. The exponents r(q) are indepen- 
dent of the ordering when L > L\ (i.e. when the ordering 
(b) becomes meaningful) and the small fluctuations from 
sample to sample are removed by ensemble averaging. 

The corresponding Renyi dimensions 

d{q) = r(q)/(q - 1) 

are shown in Fig. 4 for different L and L\, using ordering 
(a) and ensemble averaging. 

For an infinite L\ (no disorder), the eigenstates are 
plane waves of momentum k a and Q 1 ^ ^ only if k a + 
kf) — k 7 — ks = 0. This gives d(0) = 2 and d(q > 0) = 1 with 
ordering (c). The dimensions calculated with ordering (a) 



Fig. 3. Left: Spatial ordering (b). I 2 (D) for \QZ S oao \ with 
a?o in the bulk of the spectrum and for L\ (indicated by the 
arrows) = 25, 200, 2500 < L = 2500. The power law behaviors 
are obtained for a given sample if 1 < D < L\. Right: Same 
samples using spectral ordering (a). The arrows indicates the 
lower scales associated to A(Li). The power laws are valid for 
L/Li <D <L. 

are close to this limit. For a finite L\, d(q) goes from the 
clean limit (L << L\) to an Li-independent regime when 
L >> L\. In the crossover regime (L < Li) the d(q) 
depend on L\. In the limit L » L\, the d{q) (using 
orderings (a) or (b)) do not depend on L and L\. For 
< q < 3, 

d(q) w 2 - Aq 

with a slope A k> 0.135. The Li-independence of A is 
shown in the insert of Fig. 4 for L\ < L up to L\ — 600. 

A multifractal distribution has scaling behavior de- 
scribed by the /(a)-spectrum, given by the relations: 



a(q) = and f(a(qj) = a(q).q - r(q). (4) 
We obtain 

/(«(<?))« 2 - V (5) 

for q < 3, i.e. a parabolic shape f(a) = 2 — (a — 2 — 
A) 2 / (AA) around the maximum 2 + A. We have mainly 
studied the first positive moments, since we are mainly 
interested by f(a(q = 2)). Indeed, when one uses Fermi 
golden rule to calculate the interaction-induced decay of 
a non-interacting state, one needs to know the density of 
states directly coupled by the second moment (q = 2) 
of the hopping term. The fractal dimension of the sup- 
port of this density is given by f(a(q = 2)). For greater 
values of q, there are deviations around the parabolic ap- 
proximation, indicating deviations around simple lognor- 
mal distributions. From a study of the large and small 
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Fig. 4. d(q) with ao in the bulk of the spectrum using ordering 
(a) and after ensemble averaging. Filled symbols: L — 240 
and Li = oo (diamonds), 25. 10 4 (circles) and 2500 (squares). 
Open symbols: L\ — 70 and L — 960 (diamonds), 480 (circles) 
and 240 (squares). Insert: the slope A(Li) showing that d(q) is 
disorder independent for q < 3 and Li < L. 



values of \Q2? a 1 7 one can obtain d(q — > ±00). Wc find 
d(+oo) = 1.33 and d(— 00) = 3.15, giving the limits of the 
support of f(a). 

We have also checked that our results for QZ S aaa do 
not depend on the chosen ao and studied the general case 
where \a > and \(3 > are not the same. In Fig. 5, one can 
see that the QZ S „a studied for different a give the same 
curves d(q). Using energy ordering (a) and imposing an 
energy separation |e Q — e^j| > A{L\) in order to have a 
good overlap between the fixed states a and (3, we find 
also power law behaviors for I q (D). The corresponding 
dimensions d(q) are given in Fig. 5, characterized by a 
slope 

A(a + (3) « A(a = (3)/2 w 0.065. 

Therefore, the multifractal character of Ql^j is less pro- 
nounced when I a >^ \[3 >, but remains relevant. 

So far, we have discussed the hopping terms of the 
general iV-body problem. Wc now discuss how our results 
modify previous assumptions for two interacting parti- 
cles (TIP). As pointed out by Shepelyansky, the interac- 
tion induced hopping mixes nearby in energy TIP states 
\a(3 >= d+ t d^\0 >. The decay width P [8, 9,' 10] of a TIP 
state \a(3 >, built out from two one particle states local- 
ized within Li, can be estimated using Fermi golden rule. 
If one assumes RMT wave functions inside L\ for the one 

±U.L 1 3 ^ 2 couple 



particle states, (case (iii)) the Q lS 



a/3 



the TIP state \aj3 > to all the TIP states |7<5 > inside L\. 
Around the band center, they have a density pi(L\) oc L\ 



Fig. 5. d(q) calculated for L\ — 70 and L = 240 using ordering 
(a). The filled symbols correspond to a — (5 — ao, averaged 
over 10 consecutive ao chosen in various parts of the spectrum: 
down triangle (130 < a < 140), square (120 < a < 130), 
diamond (145 < a < 155) and up triangle (170 < a < 180). 
Empty symbols correspond to the case a /3. a is fixed and 
d(q) is averaged for a few /3 / a: circle (a = 100, 123 < j3 < 
133), cross (a = 80, 123 < (3 < 133) and plus (101 < a < 112 
and 123 < (5 < 133) 



and Fermi golden rule gives 



U 2 U 2 

r(E « o) cx -^^(£1) = Ti 



(6) 



. We have shown that all the TIP states which can be 
coupled by the interaction within the localization domains 
are not equally coupled. Since the square of the hopping 
terms appears in the Golden rule, our multifractal analysis 
gives a reduced effective TIP density pf oc L{ {a{q - 2 >> 
which should replace the total TIP density p 2 (Li). The 
resulting expression 



r „ rx ^!r/(«(9= 2 )) 

-t afi OC Jz L l 

can be compared to the direct numerical evaluation: 



r a (3 



U 2 Y,\Ql S e\ 2 ^c 



7<5 



e/3 



(7) 



(8) 



of the Golden rule decay. 

In Fig. 6, we show for three different sizes L how the de- 
cay rate r aa numerically calculated using Eq.(8) depends 
on Li, for a TIP state \aa > where a is taken in the bulk 
of the spectrum. From Fig. 4, one gets a (2) « 1.52 and 
t(2) w 1.69, and Eq.(4) gives f(a(q = 2)) w 1.35. For this 
value, one can see in Fig. 6 that Eq.(7) and Eq.(8) give in- 
deed the same Li-dependence. This observed L^ 1 ' 65 law 
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TIP eigenbasis for U = 0. In a third paper [16], a study 
of the dynamics of a TIP wave packet will be presented, 
showing that the center of mass exhibits anomalous diffu- 
sion between L\ and L 2 . These three studies provide con- 
sistent and complementary observations supporting our 
claim: multifractality and criticality are relevant concepts 
for a TIP system with on site interaction in one dimen- 
sion. Our results go beyond the TIP problem and show 
that oversimplified two-body random interaction matrix 
models [17, 18, 19] which ignore multifractality in the hop- 
ping cannot properly describe the many body quantum 
motion in Anderson insulators. 

We are indebted to S. N. Evangelou for very useful 
comments. 



Fig. 6. r aa {L\) for three different sizes L. Dashed line: y — 
2.55L" 1 - 65 . Doted line: y = 0.5L" 1 

clearly differs from the L^ 1 law implied by the RMT as- 
sumption (case (iii)). We can also see that r aa does not 
depend on L when L\ < L, since there are no significant 
hopping terms for range larger than L\. 

Another interesting issue is the enhancement of the lo- 
calization length L 2 , which is induced by the interaction 
and characterizes a restricted set of TIP states which have 
a sufficient overlap to be re-organized by a local interac- 
tion [11, 12]. Using the Thouless block scaling analysis [13], 
one finds jfi oc (±U / L 3 / 2 p 2 (L 1 )) 2 . If the density p 2 (L 1 ) 
of states coupled by the interaction is the total TIP den- 
sity for a size L\, one finds the original estimate [4,13] 
L 2 oc L\. The multifractality yields a reduced effective 
density p| ff oc L^ a( - q - 2 ^ instead of the total TIP density. 
Since the contribution of TIP states \a(3 > with a / (1 
dominates, we use f(a(2)) = 1.75 valid when a ^ (3 and 
we find Li oc L\ 5 . This Li-dependence is in agreement 
with recent numerical results [2,14]. So there is an en- 
hancement, though weaker than the original estimate [4] 
(L 2 oc L\), due to the multifractal distribution of the hop- 
ping terms. 

In summary, we have studied how one particle dy- 
namics (one dimensional localization) can affect the many 
body problem through non trivial properties of the dis- 
tribution of the two-body interaction. In a clean system, 
one has f(a(2)) = 1 and the density of states which are 
effectively coupled by the interaction is the one particle 
density p\ oc L. The disorder, as it is well known, en- 
hances the effect of the interaction, since the effective den- 
sity pf oc L{ [a(2)) , with 1 < f(a(q = 2) < 2 for L = L x . 
This enhancement of the density of states coupled by the 
interaction inside a system of size w L\ is nevertheless 
smaller than the one (p 2 oc L\) given by fully chaotic one 
body states inside their localization domains. In a sec- 
ond paper [15], a study of the TIP spectral fluctuations 
will be presented, showing that statistics is critical (as for 
the one body spectrum at a mobility edge) if U is large 
enough, accompanied by multifractal wavefunctions in the 
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